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Abstract. In this paper we present two concrete models of non-perfect fluid with bulk viscosity to interpret 
the observed cosmic accelerating expansion phenomena, avoiding the introduction of exotic dark energy. 
The first model we inspect has a viscosity of the form C = Co + (Ci ~ (2<l)H by taking into account of the 
decelerating parameter q, and the other model is of the form C^" = C^o + CiH + CzH^. We give out the exact 
solutions of such models and further constrain them with the latest Union2 data as well as the currently 
observed Hubble-parameter dataset (OHD), then we discuss the fate of universe evolution in these models, 
which confronts neither future singularity nor little/pseudo rip. From the resulting curves by best fittings 
we find a much more flexible evolution processing due to the presence of viscosity while being consistent 
with the observational data in the region of data fitting. With the bulk viscosity considered, a more 
realistic universe scenario is characterized comparable with the ylCDM model but without introducing the 
mysterious dark energy. 



1 Introduction 

The independent discovery respectively in 1998 and 1999 
indicates that the current universe is in accelerating ex- 
pansion [1] . To accommodate this surprisingly exotic phe- 
nomenon, dark energy, varieties of models are proposed. 
The basic idea of dark energy comes up in the context of 
supposing the general theory of relativity works precisely 
well in cosmological scale, a perfect fluid with effectively 
large enough negative pressure is required to speed the 
universe expansion up. According to the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) 7-year dataset anal- 
ysis [2] it makes up about 72.8% of the universe's con- 
tents. In the past decade, many attempts are made to 
understand the accelerating mechanism of dark energy, 
the most mysterious component of the universe hitherto 
as envisioned. The ylCDM is such a model and receives 
great attention because it is well consistent with the ob- 
servational data although suffering from some serious fun- 
damental physics problems [3] . Another popular theory of 
dark energy is referred as quintessence [4]. However, an 
alternatively instructive idea is that the general theory of 
relativity may fail in large scale, therefore modification of 
gravity theory should be introduced to drive an acceler- 
ating phase of universe expansion. The f{R) gravity, for 
example, which generalizes the Einstein-Hilbert action by 
invoking an arbitrary function of Ricci scalar, is of this 
class [5]. More details and references are available in the 
recent reviews [6]. 
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In the context of perfect fluid, models with different 
equation of state are extensively studied (see Ref. [7]), 
such as Chaplygin gas, generalized Chaplygin gas, inho- 
mogeneous equation of state, and barotropic fluid dark 
energy, etc. However, the perfect fluid assumption is as- 
sertive since it suggests no dissipation, which actually ex- 
ists widely and intuitively plays a critical role in the uni- 
verse evolution, especially in the early hot stages. To be 
more realistic, models of imperfect fluid are invoked by 
introducing viscosity into the investigation. The simplest 
theory of the kind is constant bulk viscosity, regarded as 
equivalence to constant dark energy with cold dark mat- 
ter model. A widely investigated case is that bulk viscosity 
with the form as a linear function of the Hubble parame- 
ter, which is proved to be well consistent with the observed 
late-time acceleration and can recover both the matter and 
dark energy dominant eras. Another case is that bulk vis- 
cosity behaves as a parameterized power-law form with 
respect to the matter density, which can be shown to be 
similar to the Chaplygin gas model [8]. Furthermore, vis- 
cosity term contains red-shift and higher derivatives of 
scale factor is also considered as a more general theory [9] . 
On the other hand, as subsequence of viscosity and dissi- 
pation, turbulence is considered as a cosmic component in 
Ref. [10]. In this paper we focus on the bulk viscosity and 
propose two such new models, which can be dealt with 
analytically and performed well when comparing with ob- 
servational data. 

It is well known that significant number of dark en- 
ergy models suffer from the finite-time singularity prob- 
lem. The classification of the (four) finite-time future sin- 
gularities has been proposed by S. Nojiri et al. (2005) [11], 
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while the big rip (Type I) and sudden rip (Type II) are 
discussed previously in Ref. [12] and [14], respectively. Re- 
cently, a novel concept has been proposed in Ref. [15], the 
so-called "Little Rip" , in which for some kinds of models 
the non- viscous dark energy density increases with equa- 
tion of state parameter co < —1, but that a; ^ — 1 asymp- 
totically is required to avoid the future singularity, thereby 
leading to a scenario that bound structure is disassembled 
(see Ref. [16] for more concrete models of this class). Like- 
wise, "Pseudo Rip", which is introduced as another case 
of the rips, occurs when an upper bound Finert where F 
refers to the incrtial force is confronted during the increas- 
ing of dark energy density with scale factor. Several mod- 
els have been studied as this type [17]. Further more, mod- 
els reconstructed to realize or avoid such rip/ singularity 
in scheme of modified gravity are also studied recently in 
Ref. [18]. The presence of rips causes destruction of cosmos 
structure, therefore it is an important aspect for a cosmo- 
logical viable model building that a natural scenario can 
be yielded to cure such problem [13]. In the scheme of 
viscous cosmology, singularity as a vital feature attracts 
great attention (see most of the models in Ref. [7] [8]). 
Recently, the rips confronted in viscosity models are also 
studied in Ref. [19]. Generally, it implies that the presence 
of particular viscosity could mitigate the singularness of 
cosmos quantities in some extends. 

In this letter, we propose and investigate two different 
viscosity models, both of which arc constructed under the 
assumption that our universe is bound more tightly there- 
fore the bulk viscosity varies with respect to not only the 
conventional momentum term H, but also the accelerating 
status {H, model I) and total energy {H^, model II). Ar- 
rangement for this article is as follows. In the next section, 
we derive the explicit solution of the models, respectively, 
and calculate analytically the corresponding cosmic quan- 
tities. In Section 3, wc briefly review the statistics analysis 
employed for constraining model parameters by SNe Type 
la date and OHD dataset. In Section 4 figures and tables 
are listed as our results, in which the comparison with 
the ylCDM model is also performed. In Section 5, as our 
conclusion, we discuss the result of best-fitting and the 
fate of universe evolution accordingly by means of the sin- 
gularity and rip analysis, as well as the advantages and 
disadvantages of such viscosity models. For the benefit to 
the related astrophysics and cosmology community we list 
the latest 18 OHD in the appendix. 



2 Cosmology Models of Imperfect Fluid 

In the scheme of imperfect fluid, the energy stress tensor 
reads, 



viscosity, and = -Khf^^iT^ + TWUf,,,,) the heat flux 
density four-vector with k the thermal conductivity. 

In the case of thermal equilibrium, = 0. More- 
over, the term of shear viscosity vanished when a com- 
pletely isotropic unverse is assumed. In the isotopic and 
homogeneous Friedmann-Robertson- Walker (FRW) met- 
ric, 



dr^ 



1 — fcr^ 



r^dn^ 



(2) 



where k = —1, 0, 1 being the curvature parameter, Eq. (1) 
can be rewritten as 



T^u = pU^U, + {p-eQh^,. 



(3) 



Here 9 = 3-ff where = ^ is the Hubble parameter and 
the dot denotes differential with respect to cosmic time t. 
Therefore the corresponding Friedmann equation, consid- 
ering the case of flat space-time, is of the form 



H + H 



9 47rG , 



(4) 
(5) 



as well as the equation of energy conservation for a com- 
plete dynamics system. 



p+{p + p)e = 0, 



(6) 



where p = p — 6( is the effective pressure. If we assume the 
matter presented is cold thereby pressureless, p will have 
a simple form of — ^C- So far we can principally solve the 
matter density p with respect to cosmic time t when the 
explicit form of ( is specified. 

Combining equations (4), (6) and the current value of 
Hubble parameter Hq, we obtain the following dimension- 
less equation, 

= ih, (7) 



Ho 



- 2 H„C 



the dimcnsionless 



in which h = -n- and f ^ ^ ^• 
bubble parameter and viscosity, respectively. Using the 
relation 

— = (8) 
dt a d In a 

and let ' denote the differential respect to conformal time 
In a, we obtain a dimensionless equation. 



(9) 



T^,^= pU^,U^ + {p-eC)h^.-2r]af,^ + Qf,U^ + Q^Uf„ (1) 

where p is the mass density, p the isotropic pressure, = 
(1, 0, 0, 0) the four-velocity of the cosmic fluid in comoving 
coordinates, /i^j, = -I- J/^t/^ the projection tensor, 6 = 
^il = ^'^^fi^ia-./J) = C^r^ the expansion scalar, cr^^ = — 
j/i^j/^ the shear tensor, ( the bulk viscosity, r] the shear 



2.1 Model I: Viscosity with Decelerating Parameter q 

In this subsection we consider the model by taking into ac- 
count of the decelerating parameter q. The bulk viscosity 
reads, 

C = Co + (Ci - C2q)H, (10) 
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in which q = — is referred as the decelerating param- 
eter, and Ci)C2,C3 are constant, respectively. By introduc- 
ing decelerating parameter q into the viscosity, the model 
now can take account of the role played by the transition 
of cosmic evolution phases. 

By transformation Co = ^^^o, Ci = ^ (Ci - 6) 

and C2 = ^^P-$2, we obtain the dimensionless form of 

-"0 

viscosity, 

e = & + ei/i + 6/i', (11) 

where the prime denotes Considering the condition 

/i(ao) = 1, Eq. (9) can be solved analytically when 7^ | 
and ^2 7^ 1 as 



h{a) = 



2^0 



1 



2^0 



(a/ao) . (12) 



3-2^1 V 3 - 26 
After integration we obtain the corresponding scale factor, 

(13) 

Hence the Hubble parameter with respect to t reads 



H{t) = Ho 



2^0 



3 - 2^ ^j^Hoit-to) _ 3-2go-2;i ' 



(14) 



3-2{i 



and the decelerating parameter q is 



Substitute Eq. (14), (15) to Eq. (10) we obtain the evolu- 
tion of bulk- viscosity, 

_ 127rG / 3go Co(26 ~ 36)(3 - 2^0 - 2^0 
Ho V3-2ei (3-2^)2(1-6) 



(16) 



^jil^Hoit-to) _ 3-2{o^-2Ci 



3-26 



2.2 Model II: Viscosity with H'^ 

In this subsection we consider the viscosity model with 
Hubble parameter up to quadratic order. 



C = (o + CiH + (2H^ 



(17) 



The corresponding dimensionless viscosity is given as 

^ = eo+a/i+6/i', (18) 

We should note that it is another case of stronger coupling 
between matter and dark components since the term ~ 
in proportion to velocity square (kinetic energy) takes 
into account of the energy transference. We will discuss 
the detailed feature of this model in Section 4 and 5 
The solution of this model depends on the sign of zi = 

(i-6)'-466- 



2.2.1 z\ < 

In the case of A <0, when 2^ + 4^ — 3^0 the evolution 
of universe is. 



h = 



3 - 4^ - 26 
26+46-3 



1 - 



3 - 46 - 26 
26+46-3 
1 



2;i+4{2-3 



-tan In f- 

2\/—A ^ oo 



• (19) 



Once more, integration is needed to obtain the behavior 
of scale factor with respect to t. However, it can not be 
shown explicitly due to the complexity. Therefore we will 
directly assign the best-fit values for the coefficients and 
calculate the corresponding cosmic quantities numerically. 



2.2.2 Z\ > 

In the case of /4 > 0, an alternative expression is yielded 
when 2/i ^ 26 +46 + 3 and 6 7^ (iJ^ term appears) 

as 



(3-26)-27Z 7Z 2VZ + (26 + 46 + 3) 
6 2aM-(26+46 + 3) 




2y/A + (26 + 46 + 3) 
2v^-(26+46+3) 

for simplicity we rewrite this formula as 

h 



. (20) 



1 



66 , 66 

6 - 1 6 - (a/ao)«= 



(21) 



by redefining 6 = ^^+^^^^+^^^+3) and 6 = From 

2v^-(2Ci+452+3) C2 

^ = hHo an expression of a{t) is obtained by integration, 
which is 



a{tf 



1-«1«2-1 



= Ci(6-66-i) 



(6 - a{t)(^ 

'(6-66-i)(6-6-i) 



X exp 



6-1 



Ho{t-to)\, (22) 



where Ci is an integrating constant determined by the 
condition a(to) = oo- We cannot solve a{t) explicitly due 
to the complex form. Nevertheless, we can study the fu- 
ture behavior by appropriate approximation according to 
the numerically fitting result. When a{t) grows with time 
flying so that |a(t)^^| 6) Eq. (22) can be easily calcu- 
lated as 



I,. / (6-66-i)(6-6-i) ^ , ' 

a{t) ~ao exp — ^ — ^ j Ho{t - to) 



(6-i)(6+6-266) 



(23) 



4 



Xin-he Meng, Zhi-yuan Ma: Rip/singularity free cosmology models with bulk viscosity 



Under the same approximation the Hubble parameter with 
respect to t is constant, 



H{t) —Hiate-time 



(a -66-1) (6 -6-1) 

(6-i)(6+6-2a6) 



Ho 



(24) 



Thus /io is automatically minimized as /Uq = ^ by calcu- 
lating the following transformed x^: 



xIj,M = Ail!}) - 



C ' 



(32) 



3 Data Analysis 

3.1 Type la Supernovae 

With its cixtraordinary property of uniform absolute mag- 
nitude, observation of Type la Supernovae suggests a way 
for inspecting the history of universe by constructing the 
relation between the red-shift z and luminosity distance. It 
is believed that the SNe Type la observation can provide 
the most direct evidence of cosmic accelerating expansion. 
In this paper we use the Union2 SNe la dataset [20] for 
best-fitting, which compiles 557 SNe la covering the red- 
shift range z = [0.015, 1.4], as shown in Fig. 1(a). It ex- 
tends the Union dataset by including new data points at 
low and intermediate redshifts discovered by the CfA3 and 
SDSS-II Supernova Search projects, as well as 6 new SNe 
discovered by the Hubble Space Telescope at high z. 

To perform the statistics analysis, the theoretical 
distance modulus is defined as 



Hth = Slogio-Di, iz)+no, 



(25) 



where = Hod(^z) is the dimensionless luminosity and 
= (1 -h z) du {z), with du {z) denoting the co-moving 
distance, 

Therefore, the corresponding Xsjve function is calculated 
from 



XSJVe 



E 



O'obsiZi) 



(27) 




0.0 0.2 0.4 0.6 



0.8 



1.0 
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1.4 



(a) The SNe la dataset contains 557 points range 
from redshift 0.015 to 1.4. 
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(b) The OHD dataset constains 18 points range 
from redshift 0.09 to 1.75. 

Fig. 1. The observed SNe la and OHD data (with error bar) 
that we employ for best-fitting approach. 



which should be minimized consequently with properly 
choosing of /xq and model parameters d. The minimization 
with respect to /uq can be made trivially (Ref. [21]) by 
expanding xsNe as 



XsNe =A- 2noB + ulC 



where 

n 



i=l 



IJ-obsizj) - jJ-thjzj-.'d-.fio = 0) 

<^obsiZi) 

B{'d) = ^ l^obsi zj) - IJ,th{zi; I?; /xq = 0) 



i=l 



(^Ibsi^i) 



n ^ 



(28) 

(29) 
(30) 
(31) 



3.2 Observational Hubble Parameter Data 

Recently, the direct measurement of H{z) arouses much 
attention. In this paper we use the OHD dataset of 18 mea- 
surement points collected in Ref. [22], which combines the 
observational constraints on the Hubble parameter avail- 
able so far. It includes 8 points of the Hubble parameter 
in the range 0.1 < z < 1.75 from the relative dating of 
32 passively evolving galaxies [23], 2 more observations 
of the Hubble parameter at z ^ 0.5 and z ^ 0.9 ob- 
tained by high-quality spectra of red-envelope galaxies in 
24 galaxy clusters [24] , and 8 new high-accuracy estimates 
of H(z) provided by Ref. [25] for helpful compiling, ex- 
tracted from several spectroscopic surveys among a large 
sample; of 11324 early type galaxies. All the values arc re- 
ported in Table 2 in the appendix and shown in Fig. 1(b). 
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The for Observational Hubble Data is 



XoHD — 5Z 



(33) 



where h{z) = H{z)/Hq is the dimensionless Hubble pa- 
rameter. Following the same approach as XsATei the mini- 
mization with respect to Hq can be satisfied automatically 
by introducing the transformed function: 



XX)HD 



+ C, 



where 



1=1 * 



■J^ h{zi)Hobs{zi) 

^ = 2^ ^5 ^ 



i=l 



(34) 

(35) 
(36) 
(37) 



4 Results 



In this section wc perform best-fitting to the models men- 
tioned in Section 2 and plot the evolution of relative quan- 
tities, respectively. For computing feasibility, we shall rewrite 
Eq. (12), (19), and (20) with respect to redshift z and in- 
troduce alternative parameters for smoothness of the gra- 
dient calculated in minimizing x^- 



• Model I, bulk viscosity with q, 

/l(z) = (l-Al) + Al(l + ^)^^ 



(38) 



with Ai = 1 



3^ and A2 



Model H.a, bulk viscosity with iJ^, Z\ < 0, 

h(z) = Xi + (1- Xi)- ^ 7T— t-t; TT: (39) 

^' ^ a - A2 tan(A3 ln(l-t- z)) ' ^ ' 

with A, = A. = and A3 = 

Model H.b, bulk viscosity with H^, A > 0, 



h{z) = 1 



+ 



A1A2 



A1A2 

Ai + 1 ' {l + z)-^^+Xi 



(40) 




(a) Model I. A transition from negative to positive 
viscosity is shown in this figure. The negative vis- 
cosity corresponds to the era of inflation which has 

sufficiently largo q. Then the value turns positive 
and increases until approaches to a constant value 
in the future. 




(b) Model II with A < O.lt increases from negative 
to positive and ends up with zero after a summit 
at late-time. 



\ 
V\ 
V\ 
\\ 










SNe 

OHD 




\\ 
\\ 
\\ 








\ 









1.0 1.5 2.0 



-1.0 -0.5 0.0 0.5 

t-to 

(c) Model II with Zi > 0. It begins at a positive 
value then drops below zero and recovers asymp- 
totically to a lower positive value. The negative 
viscosity corresponds to the deceleration of mater 
dominant phase. 

Fig. 2. The 5 — t relation diagram of Model I, Il.a and Il.b 
when best-fit values of parameters are given. 



with the ^ in Eq. (21) replaced by A. 

The results of data fitting is listed in Table 1, while the 

calculated evolution of cosmic quantities plotted respec- 
tively in Fig. 3, 5 and 6. For comparison, we also plot the 
evolution of ^CDM model, H"^ = f2mO'~^ + 1— Om, with 

Qrn = 0.27. 



5 Conclusions and Discussions 

In this paper we continue our previous work on bulk vis- 
cous (dark fluid) cosmology by focusing on two explicit vis- 
cous models and studying the viability of accommodating 
the observed accelerating expansion of the universe with- 
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2.0 



-0.5 0.0 0.5 

t-to 

(a) The a — t relation diagram. Basically there are 
three phases during the evolution: An exponen- 
tially inflationary scenario at the very beginning, 
followed by a decelerating phase (matter dominant 
era), eventually enters into the accelerating expan- 
sion (dark energy dominant era). Wc can sec the 
results correspond to TlCDM model very well and 
constraint by SNe performs better than OHD. 
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0.5 
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1.0 1.5 2.0 



(b) The H — t relation diagram. The curves of SNe 
and ylCDM still stay close while the result of OHD 
deviates slightly. 




-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

t-to 

(c) The q — t relation diagram. In this panel we see 
significant distinction between the proposed model 
and ylCDM in the early time of universe. Unlike 
the initial value 0.5 in ylCDM, the decelerating pa- 
rameter of this viscous model starts at a relatively 
large value and rapidly falls to and crosses over 
zero, followed by a smooth evolution to -1. 

Fig. 3. Resulting evolution of cosmic quantities of Model I 
when best-fit values assigned. Result of ACDM model is also 
plotted for comparison. 





SNc 







10 
t-to 



15 



20 



Fig. 4. The far-future a — t relation diagram of Model II with 
A <Q. The scale factor in this case behaves exotically that it 
climbs up to an upper bound at time tc, then remains constant. 
We will discuss the feature in Section 5. 
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11. a < 


11, A X) 




Ai = 


0.15312 


Ai = 0.11577 


Ai = 0.039259 


SNe 


A2 = 


= 2.5566 


A2 = 0.76136 


A2 = 3.2000 






= 542.38 


A3 = 0.88993 
xLn = 542.15 


xLn = 542.59 




Ai = 


0.34146 




Ai = 0.041795 


OHD 


Aa = 


-- 1.7833 


N/A 


A2 = 3.2641 




Xmin 


= 12.513 




xLir. = 12.005 



Table 1. Best-fit results of Model I, Il.a, and Il.b 



out introducing the mysterious dark energy. We employ 
the Union2 data and OHD dataset to generate the best- 
fit values for the models, and then accordingly draw the 
corresponding evolution diagrams, respectively. As com- 
parison, wc also demonstrate that the resulting evolutions 
of such models are precisely consistent with /ICDM model 
in the region of observational data. Both of our models 
are proposed in the context of assumed tight coupling be- 
tween the cosmic components, and by arousing a time- 
varied bulk viscosity as constructed, we see the cosmic 
evolutions of such models can be greatly different from 
that suggested by the ylCDM. 

Wc present our main conclusions and discussions in 
order as following. Model I, C = C(<?) -f^)) mimics the uni- 
verse evolution of ylCDM model perfectly. We can see from 
Fig. 3(a) that the unverse has an exponentially inflation 
after the Big Bang, and soon enters into deceleration when 
in matter dominant era, followed by a late-time acceler- 
ating expansion. The monotonic increasing bulk viscosity 
(Fig. 2(a)) with respect to both q and H provides a con- 
nection among the pressurelcss matter therefore behaves 
as effective dark energy in the late-time evolution. We 
should notice that although the age of universe is smaller 
in this model than that of ylCDM, it shares the most sim- 
ilarities with those of the TlCDM model, which makes this 
kind of bulk viscosity component model a successful sub- 
stitution for dark energy. 

Model II, C = C{H,H'^) as another case of viscous 
cosmology model, suggests two alternative evolutions due 
to different values of parameters. In the Case of Z\ < 0, 
we obtain a bound universe that the expansion ends at far 
future t = tc, when H and H both vanishes (static then) 
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3.5 
3.0 
2.5 
^ 2.0 
1.5 
1.0 
0.5 



t-to 

(a) The a — t relation diagram. We can see in this 
case the evolution of scale factor coordinates with 
ACDM model in the range of observation data, and 
the evolution curve has a similar trend with Model 
I (Big bang, inflation, deceleration, and accelera- 
tion). Nevertheless, there arc two main differences: 
The age of universe is much shorter than that pre- 
dicted in ylCDM; The future behavior is widely di- 
vergent (see Fig. 4). 




t-to 

(b) The H — t relation diagram. Instead of a mono- 
tonic descending followed by an asymptotical con- 
stant value as suggested by ylCDM model, the evo- 
lution of Hubble parameter in model decreases 
more dramatically at early time, and goes towards 
zero in the future. 



1.5 

1.0 

0.5 

0.0 
-0.5 
-l.OL 




-0.5 0.0 



0.5 
t-to 



1.0 



1.5 



2.0 



(c) The q — t relation diagram. In accord with the 
future halt of expansion, the run of q is essentially 
different from that of ACDM. 

Fig. 5. Resulting evolution of cosmic quantities of Model II 

with best-fit values assigned when Zi < 0. Wc only plot that 
constrained by SNe data, because a poor confidence of esti- 
mated parameter values is confronted when OHD dataset is 
employed. 



3.5 
3.0 
2.5 
= 2.0 
^ 1.5 
1.0 
0.5 
0.0 



OlID 
LCDM 



-1 



t-to 



(a) The a — t relation diagram. We find no big 
bang in this case, instead of which the evolution 
of scale factor performs as two de Sitter-like ex- 
pansion combined via a transition phase. However, 
in region of observational data, the model corre- 
sponds to ylCDM very well. 




SNe 
OHD 
LCDM 



t-to 

(b) The H — t relation diagram. We can see clearly 
the three phases of the evolution in this panel, be- 
ginning with a relatively larger value of H, then 
going down smoothly during the transitionary era, 
and ended up with a smaller value of H . 




0.0 0.5 1.0 1.5 2.0 

t-to 

(c) The q — t relation diagram. In accord with the 
three phases, the decelerating parameter q goes 
from negative to positive in the pre-big bang era, 
then runs very closely to that of ylCDM model in 
the data region and future. 

Fig. 6. Resulting evolution of cosmic quantities of Model II 

when Z\ > 0. In this case wc obtain a universe that has no big 
bang. It acts as de Sitter expansion with different values of H 
in both the very past and far future. The range of data fitting- 
locates in the transitionary era between the two phases, which 
also mimics precisely the behavior of ylCDM in that region 
when the best-fit values are assigned. 
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(Fig. 4). However, it fits the data quite well in correspond- 
ing data region as seen in Fig. 5(a). The halt confronted 
in the far future is the result of the viscosity with H'^ that 
drives (7) to a trivial solution. From the numeric calcula- 
tion we estimate the value of that time, tc — to ~ 70.2Gyr. 
In the Case of Zi > 0, a no-big-bang unverse is obtained. 
It performs as a dc-Sitter expansion in both the very past 
and very future with different values of H. The range of 
data fitting locates in the transitionary phases between the 
two eras. In the region it mimics the observed acceleration 
successfully. From Fig. 6(b) and Fig. 2(c) we see that the 
values of H, q and ^ go to approximated constants when 
tc-toc^ 19.4Gyr. 

We are able to investigate the universe fate of our cos- 
mology models by means of the rip or singularity. Singu- 
larity behavior is classified in Ref. [11] as follows: 

• Type I ("Big Rip"): For t -> ig, a -> oo, p — >■ oo, and 

\p\ oo; 

• Type II ("Sudden"): For i — >• t^, a — >■ a^, p — )■ ps, and 
\p\ oo; 

• Type III: For t ^ tg, a a^, p — > oo, and \p\ — > oo; 

• Type IV: For t ^ tg, a ^ Us, p ^ 0,and \p\ —J- 0, and 
higher derivatives of H diverge. 

In our cases, an effective pressure is arisen as a result of 
the presence of bulk viscosity, reads p = —3H(. Hence 
the evolution of \p\ is similar to that of bulk viscosity in 
the future because of the asymptotic constant for Hub- 
ble parameter. According to Fig. 2, 3, 5 and 6, none of 
the finite singularity is confronted, therefore viscosity in 
our model as a practical substitution of dark energy can- 
didate is better than some other theories (e.g. phantom 
dark energy). 

The Rip behavior is another aspect we interest in. 
Three types of Rips are discussed in Ref. [17] according 

to the evolution of inertia! force on a mass m as seen by 
a gravitational source separated by a co-moving distance 
I, defined as Finert = ml{H + H'^). 

• "Big Rip": For t -> ts, Finert ^ oo; 

• "Little Rip": For t — >• oo, Finert oo; 

• "Pseudo Rip": For f — )■ oo, Finert never goes infinite 
but has an upper bound. 

Rips occur when the bound structures are disintegrated 
by the cosmic expansion, thereby being an important fea- 
ture in testifying models. In our models, the incrtial force 
is calculated by Finert oc H'^ + 9H(, which never varies 
towards infinite during the evolution, i.e., they are free 
of little rip. What's more, although in model I and II. a 
there are maximums of Finert, they are merely fail pseudo 
rips because of the small value that cannot disassemble 
the bound structure of the universe, i.e. they are free of 
pseudo rip. 

Our real universe is filled with viscosity media, there- 
fore cosmology based on imperfect fluid is proposed rea- 
sonably, which is proved to be viable in interpreting the 
observed cosmic accelerating expansion. We see that for 
the theories we study in this paper, the presence of bulk 
viscosity allows for a much more flexible evolution desti- 
nation, when avoiding any type of singularity or rip. On 



the other hand, in this work we have only pressureless 
(cold dark) matter and an effective pressure provided by 
viscosity term, therefore it is hard to recognize Om from 
the total matter density, as they arc coupled tightly, if 
the form of viscosity or its evolution behavior is complex. 
This explains why we discard to use CMB-shift and BAO 
data for joint constraint analysis. We also show that, as 
a promising probe for cosmological fitting, observational 
Hubble data performs very well and gives reliable results. 
It is believed that with much more high-precision data 
points released or published, the OHD dataset will play 
an even important role as an independent and direct mea- 
surement. 

Instead of adding bulk viscosity to T^y, the R.H.S. 
of Einstein equation, modification on the gravity theory, 
the left side of Einstein equation or the gravity sector, is 
another approach of accommodating the cosmic late-time 
acceleration expansion. It is possible to be reconstructed 
effectively from the left side to the right side of the equa- 
tion, which is also intriguing and we will be present that 
explorations elsewhere. It is also expected in the future 
work to give additional constraints on the parameters of 
viscosity models from cosmic fluid dynamics. 
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Appendix 

The list of 18 OHD data used in this paper. 
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